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►> Abstract: In this paper we compare the prediction for deep inelastic scattering from N=4 SYM with the 
^O HERA experimental data. The paper conveys two results. The first is the message that N=4 SYM is able 
£C) to describe the DIS data with very good accuracy (x 2 /d.o.f. < 1.5) in the region of Q 2 = 0.85 -r 60GeV 2 
^ with 2/vA = 0.7-^0.8 . The second is that the value of string coupling constant g s turns out to be so small 
t ■ that none of saturation effects will be visible in the region of accessible energies, including the maximal 
§ energy of the LHC (W = 14 TeV). 
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1. Introduction 

It is well known that N=4 SYM together with AdS/CFT correspondence allows us to study theoretically 
the regime of the strong coupling constant jjj . For the first time we have a theory which leads to the 
main ingredients of the high energy phenomenology such as the Pomeron and the Reggeons, in the limit 
of strong coupling. On the other hand, N=4 SYM with small coupling leads to normal QCD like physics 
(see Refs. y, |[) with OPE and linear equations for DIS as well as the BFKL equation for the high energy 
amplitude. 

The Pomeron which appears in N=4 SYMfJ has the intercept and the slope of the trajectory that are 
equal to 

ajp(O) = 2 - ^v = 2 " P\ a 'p(°) = °- ( L1 ) 

in the limit of p <C 1. First, we would like to recall that N=4 SYM has a simple solution for the following 
set of couplings: 
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a YM = — — -; R = a'2 \i; g s < 1; but A > 1 (1.2) 
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where R is the radius in AdS$- metric: 

ds 2 = \ I dz 2 + J2 dx i ) = ^2 ( dz2 + dxu.dx") (1.3) 

with fi = 0,1,2,3. 

One can see that at large A the Pomeron intercept is close to 2 and, therefore, the exchange of 
the Pomeron gives almost real amplitude. Indeed, the unitarity constraint in this case looks as follows 
|,|,|,0,1,1,0 

Im A (s,b;z,z) = \A (s, b; z, z') \ 2 + O ( p = -^= J (1.4) 



Eq. ( |l.4|) means that the contribution of the multiparticle production is small for the strong coupling and 
main source of the total cross section is originated by elastic and quasi-elastic ( diffractive) processes when 
the target (proton) remains intact. Such a picture not only contradicts the QCD expectations)!!, 12, [D|, 



14, 15, |16|] , but also contradicts available experimental data. 

On the other hand, the main success of N=4 SYM has been achieved in the description of the multiparti- 
cle system such as quark-qluon plasma and/or the multiparticle system at fixed temperature [17, 18, 191 2(1]. 



Therefore, we have either to find a new mechanism for multiparticle production in N=4 SYM (see an 
attempt in Ref . pi] ) or to assume that A is not very large (say p = 0.5 -f- 0.8). It should be noticed that 
even at p = 0.8 A is rather large « 6. 

The first goal of this paper is to find the range of A that can describe the deep inelastic(DIS) data 
from HERA. We believe that correction to Eq. ( |1.1|) is proportional to 1/A and A ~ 6 could lead to a 
good in the description of the experimental data. It can give a sizable cross section for the multiparticle 
production. Such an approach will be a N=4 SYM motivated model which will be able to provide a guide 
for a theoretical approach to QCD in the region of strong coupling. 

It was shown in Refs.H, 1 0, ^, ^, [1Q] that we face the saturation phenomena in N=4 SYM at low 
x for DIS. The physics of the saturation looks very similar to the saturation phenomena in high density 
QCD with one essential difference: the saturation in N=4 SYM we can theoretically describe in very simple 
fashion based on the eikonal formulae. This approach can be easily generalized for scattering both dense 
and diluted systems. Using these formulae we can learn what can happen in the region of very low photon 
virtualities and find parameters such as g s that govern the strong interactions. 

The second goal of this paper is to find parameters of N = 4 SYM that characterize the strength of 
interaction in the large coupling limit from the comparison with DIS data and derive the estimates for the 
expected saturation effects. 

N=4 SYM being conformal invariant theory has only massless particles and leads to the scattering 
amplitudes that fall as a power of b at large values of impact parameters (b). This decrease results in the 
power- like dependance of typical impact parameters in the amplitude. In particular, for hadron-hadron 
scattering these typical b oc s 1 ' 3 [10, 22| which contradicts the Froissart theorem [23|. We have to go 



beyond of N=4 SYM and discuss the string theory, conformal limit of which is N=4 SYM, to restore the 
logarithmic behaviour ( b oc Ins) of the typical impact parameters. The third goal of this paper is to find 
out how the b behaviour influences the description of the experimental data. 

In the paper we compare the N=4 SYM formula for deep inelastic structure function F 2 which we derive 
in the next section, with the HERA data for low x region (x < 0.01). Since the physical meaning of the fifth 
coordinate z (see Eq. ( |1.3|) ) is the typical size of the colliding particles, DIS gives an unique opportunity 
to check the predicted behaviour on z. On the other hand, it is known that the energy dependance of DIS 
is rich and, in particular, F 2 oc x with A = 0.1 -j- 0.5 for Q 2 = 0.1 4- 27GeV 2 , respectively. Therefore, 
we have the set of the experimental data both for checking energy and z dependance of the scattering 
amplitude, especially because the experimental errors are so small that it is a challenge to describe the 
data in any theoretical approach (see section 3 of this paper). 

The paper conveys two results. The first is the message that N=4 SYM is able to describe the DIS 
data with very good accuracy (x 2 /d.o.f. < 1.5) in the region of Q 2 = 0.85 -=- 60 GeV 2 with p = 0.7 -r 0.8 ( 
see section 3). The second is that the value of g s turns out to be so small that none of saturation effects 
will be visible in the region of accessible energies including the maximal energy of the LHC (W = 14 TeV). 
DIS data can be described both in conformal N=4 SYM and taking into account non-conformal corrections 
in b dependance. 

However for description of proton-proton scattering we need corrected b dependance (see section 3). 



The main result of Ref. [22 1 that it should be the other source of the multiparticle production than N=4 



SYM remains even for p = 0.7 -f- 0.8. 

2. High energy Scattering in N=4 SYM 

2.1 Pomeron exchange 

As has been mentioned there exists the Pomeron in N=4 SYM with the parameters of its trajectory given 
by Eq. (|1 . 1|) . The exchange of this Pomeron leads to the following contribution to the scattering amplitude 
(see Fig. ||-a): 

_ g 2 2 {z lZ2 s) l - p ln (l + « + V u ( 2 + «)) / In 2 (l + n + ^/u{2 + uf 

= ~t — 1 — ~t~ ^ r — 1 , Gxp I — - 

47rUp J V"(2 + u) ^ p ~ ln 3 (ziZ2s) y P ln( Zl z 2 s) 

where 

u = and b is the impact parameter in the scattering amplitude (2.2) 

2z\z 2 

One can see that Eq. ( |2.1| ) is very similar to the expression for the exchange of the BFKL Pomeron fl4|1 
in which the sizes of the interacting dipoles are replaced by z\ and z 2 and in which Ajp = 2 — ajp(0) and 
the diffusion coefficient are equal. It should be recalled that Eq. ( |2.1| ) describes high energy scattering in 
the kinematic region where z\z 2 s S> A 3> 1. 




2.2 Eikonal formula 

Since the Pomeron intercept is larger than 

1, considering the high energy scattering 

z 2 we cannot restrict ourselves by the exchange 

„. m , of one Pomeron. It is well known that in 

Fig. m-b 

N=4 at small coupling as well as in per- 

Figure 1: It is shown the one Pomeron ( reggeized graviton) ex- turbative QCD the problem to take into 
change in Fig. -a and the eikonal rescattering (Fig. @-b) for N=4 account all Pomer0 n exchanges and the Pomeron 

interaction is a very difficult problem that 
has been only partly solved in high den- 
sity QCD (see Refs. 0, [jj, [0| @, |l|, ||). However, for N=4 SYM with large coupling the situa- 
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tion turns out to be much simple and at small values of p the amplitude can be found in the eikonal 
approximation |], [7], |8], [H| (see Fig. |lj-b) , namely, 

A(s,b;zi,z 2 ) = ill - exp (iA(s,b;z 1 ,z 2 )) j 

The total cross section that we are going to discuss is proportional to imaginary part of the amplitude and 
can be written in the form for virtual photon- proton scattering in the form 

(rtot{l*+p) = (2-3) 

2 / d 2 b / dzidz 2 ® r (zi) ® woton (z 2 ) |l - cos nV c 2 Rei (s, b; z\, z 2 ) J exp (-N^ImA (s, b; zi, z 2 )j } 

Functions $ 7 * (z±) and ^proton ( z %) describe the probability for virtual photon and proton to have the size 
z\ and z 2 , respectively, and we will discussed them in the next section. 

2.3 * 7 * and * pro ton 

At low x the DIS on the boundary can be expressed through the dipole -proton cross section 

o~tot{l*P) = I d 2 r± P 7 * (r±) a tot (dipole-proton;r_L; x) (2.4) 



where the probability to find a dipole with the size r± P 7 * (rj_) is equal to [24] 

N f 

"s ^ 
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Pr (n) : = ^r E z % 2 + 0- - O 2 ] Q 2 K l {Qr±) (^) 

1 



where £ is the fraction of the energy that is carried by the quark, Zf is the fraction of the electric 
charge for the quark of flavour /, a em is the electromagnetic fine constant. 



To find <£ 7 * we need to generalize the wave function of the photon (Kq (Qr±)) on the boundary to the 
wave function in the bulk. We can reconstruct this wave function using the Witten formula p5|, namely, 



*-,. (r,i) = (2.6) 

^rry / dV { ^TW^f T 9 * W) with A± = 5 (" ± v ^^ 1 ^ 

where \& (rj_) is the wave function of the dipole inside the photon on the boundary. Using Eq. ( |2.6[) we can 
find $ T *(z) as 



, OtemNr. 
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iv / 1 

X>// rfCK 2 + (1 - C) 2 ] Q 2 _/ rf 2 r|^ T » (r, ^) | 2 = (2.7) 

l N f 
J dCj2zj[C 2 + (l-C) 2 ]Q d 2 rd 2 r>d 2 r» 



2vr 2 7o Z ^ ; •'/J^y ^^(A-l 

x ( , 2 + ( f- r -, )2 ) ^ (<*" J-) ( , 2 + (f- r -r,) 2 ) ^ (0rl) 

Using the formulae 3.198, 6.532(4), 6.565(4) and 6.566(2) from the Gradstein and Ryzhik Tables, 
Ref. [26 1 we can rewrite Eq. Q2.7J) introducing the Feynman parameter (£) and taking the integral over r 



and the angle between f' and r". It has the form 

JV 



/ 1 1 

Y.Z) J d([( 2 + (l-0 2 ]Q 2 Jdr 2 dr> 2 J d^K 1 (Qr' x )K 1 (Qr 1 ) 

2 (f (1 - £)(r> 2 + r- 2 ) + z 2 ) 2 + 4r'r-,e 2 (l - £) 2 



(e(i - o(t- - r ") 2 + z2 ) 5/2 (c(i - 0(r- + h 2 + -2 2 ) 5/2 

In Eq. fl2.S| ) we used that for photon m and d in Eq. (|2.6|) are equal to and 2, respectively. 



For ^proton we use the expression that has been suggested in Ref. pffi ], namely, 



^proton l z J 



/ d 2 r\\\^{ n ,z)\ 2 (2.9) 

•^ i=l 

In Eq. Q2.9| ) we assumed that a proton consists of N c colourless dipoles and each dipole interacts with other 
dipoles without correlation. We use Eq. ( |2.6j ) to find out function \£ (r{,z) . In this equation \& (r') is the 
wave function of the dipole inside the proton on the boundary. For simplicity and to make all calculations 
more transparent, we choose \P (r ; ) = Kq (qr'). The value of the parameter Q can be found from the value 
of the electromagnetic radius of the proton (q ~ 0.35 GeV^ 1 ). 



Substituting Eq. 



in Eq. (|2.9|) one obtains 



g 0-5 

^ 0.45 

()A 

0.35 
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0.25 
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0.1 
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f d 2 v d 2 r> d 2 r K Q (qr' x ) K (qr'l) (2.10) 



z 2 + (r — r') 2 / \z 2 + (r 

where AA is the norm of the dipole wave function on 
the boundary ( M = 7r/g 2 for .Ko (<Z?"±)) 

Using that 
Jo (fcr) r dr 1 



,i-A , A-l 



Figure 2: The slope dhxF^/dhxil/x) versus Q . Data 
are taken from Ref. p3|. In this paper the slope was ex- 
tracted from the data of Ref. p9| . The curves are the fit to we obtain 
the data, based on perturbative QCD, given in Ref. pa. 



( z 2 + r 2)A " r( A ) 
rdr Jq (kr) Kq (qr) = 



k 2 + q 2 



k*- 1 K^ A (kr) ; 



(2.11) 



/>oo v-2 (+\ J.2A-1 i f 

^proton (*) = 2 5 - 2A q 2 z 6 ~ 2A / ^^^ - 

Photon v ) J o (t2 + ^ 2)2 

(2.12) 



2.4 Physical observables 

We deal with the DID structure F 2 which can be written in the form 



F2(Q 2 ;x) 



Q 2 
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atoti^+p-Eq. (U)) 



For proton-proton collision we will use the total cross section written as 

otot (p + p) = / tfaidza^proton («i) ^proton {&) 
2 / (i 2 6 |l - cos(N 2 ReA(s,b;zi,z 2 )\ exp (- N 2 ImA (s , b; zi , z 2 )) } 



(2.13) 



(2.14) 



with ^proton ( z ) given by Eq. ( 2.12 ). However Eq. ( pTl] ) describes the scattering amplitude only in the 
region of high energies. For DIS we select onlt data with x < 0.02 and we use for fitting the following 
expression: 



F 2 {Q 2 ;. 



F 2 (Q 2 ;x;Eq. O) + F 2 m (Q 2 ) 



(2.15) 



where F™ (Q 2 ) is a DIS structure function at xq = 0.02 and it was considered as a fitting parameter at 
any value of Q. 



For proton-proton interaction we described the data at W = \/s > 20 GeV and add a constant a$ to 
Eq. O. 

Eq. ( |2.3| ) which is written in N=4 SYM, has power - like decrease at large valies of the impact parameter 
(b). On the other hand for not very small p the conformal symmetry of N=4 SYM is broken and we need 
to consider the string theory for which N=4 SYM is a conformal limit at small p. In the string theory the 
hadron spectrum has the lightest hadron with the mass r7 Vi ue bail = Vp/ ' wm ch leads to the exponential 



b^>m 



falldown of the amplitude A (s, b; Zi, z 2 ) at large b: A (s, b; z\, z 2 ) 
rewrite Eq. fl2.3|) in the form 



glueball] 



-> exp 



-m 



glueball ' 



We 



A(s,b; Zl .z 2 ) = A{s,b; Zl ,z 2 \;Eq. Q) e "glueball ( 2 .16) 

to take into account the mass spectrum of the string theory. The final answer for the amplitude in this 



case is Eq. (2.3) in which A is replaced by A, (A — > A). 

In this paper we check also our description of the DIS data with the experimental data on the total 
cross section for proton-proton scattering. For this observable we use the following expression: 

a proton - proton 0) = °proton-proton ( E 9- (EH)) + °"o (a) with a (s) = a i + -j= (2.17) 

The contribution ex 1/y/s corresponds to the contribution of the secondary Regge poles. <to (s) is related to 
the mechanism of the strong interaction that cannot be described by N=4 SYM or to unknown corrections 
to this theory oc 1/A. 



3. Comparison with the experimental data 



Solution 


P 


9 


A 


X 2 /d.o.f. 


I 


0.701 ± 0.004 


0.004 ± 0.001 


2 


1.42 


II 


0.75 ± 0.007 


4.019 ± 0.061 


2 


1.22 



Table 1: Fitting parameters for solution I (see Eq. (2.1)) and solu- 
tion II (see Eq. (|2.16[) ). 



Using the formulae of the previous 
section we compare the value of F 2 [Q 2 ; x) 
with the HERA experimental data in the 
region of low x ( x < 0.02) . As has been 
mentioned our main goal is to obtain two 
parameters of the N=4 SYM: p and g s . 
For each chosen value of Q 2 we introduce 
one more phenomenological parameter: the value of F 2 (Q 2 ). It should be mentioned that the value of A 
in Eq. (^Jj) for the proton as well as the value of q have to be found from the fit, but we have to recall 
that the value of q characterizes the typical scale of the non-perturbative wave function of the proton and 
can be extracted from the electromagnetic radius of the proton. 

As we have mentioned the main qualitative experimental observation is that F 2 oc (l/x) A ^ ) and 
the power X(Q 2 ) depends on Q 2 changing from A w 0.1 -f- 0.2 at low Q 2 < lGeV 2 to A « 0.4 at high 
Q 2 > 50 GeV 2 (see Fig. I) . 



At first sight we cannot reproduce such a behaviour since the exchange of the Pomeron generates only 
one power A = 1 — p. The only way out is to include the Pomeron re-scattering which could lead to the 
amplitude with the effective power that depends on Q. To our surprise we fitted the HERA data with 
small value of g which turns out to be so small that in the accessible region of energies including the LHC 
highest energy the Pomeron re-scattering does not contribute and only the exchange of the one Pomeron 
determines the amplitude. The restoration of the unitarity constraint will occur at ultra high energy. 



The quality of description one can 
see from Fig. || where the new HERA 
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data[27] are plotted at eight values of Q 2 
as function of x. The message is clear: 
the z and Ins dependence in Eq. ( |2.1| ) 
reproduces the change in A as function 
of Q shown in Fig. ||. 

The fit to the experimental data was 
made in two differrent cases: the first 
one ( solution I) corresponds to Eq. ( plj| ) 
with the amplitude A determined by Eq. ( pO| 
and the second takes into account Eq. ( 2,16| ) 
which restricts the integral over b (so- -| 
lution II). The parameters that we ob- 
tain are listed in Table 1. Solution I 
gives small value of g = N 2 g s = 0.04 
while solution II leads to rather large 
g = N 2 g s = 4. They have very close 
X 2 /d.o.f. ~ 1.5 and describe the data 
equally well (see Fig. Q). However both 
solutions cannot describe the data at lower Q (see Fig. ^|-c and Fig. ||-d). From the fit we also determine 
the function Ff 1 in Eq. ( fTl5| ) ( see Fig. |). 

It turns out that for both solutions the amplitude A is small ( see Fig. || that illustrated this fact). 
In this figure one sees that the amplitude A for the solution II reaches the value of about 0.5 but in spite 
the fact that this value does not look very small Eq. fl2.3| ) gives the value of the amplitude which is very 
close to IrmA. For the solution I the value of g is so small that it leads to a very small A. Therefore, the 
lesson which we obtain from this estimates is very simple: the data for DIS can be described in N=4 SYM 
but the non-linear (shadowing) corrections tuns out to be very small. In other words the saturation effects 
which are in N=4 SYM are very similar to the one in QCD[Q, ||, |(| [7|, [8|, [T0| , will be sizeable only at ultra 
high energy, higher than the LHC maximum energy (W = lATeV). 

These two solutions we check against the experimental data for the total cross section of proton-proton 
interactions. The comparison with the experimental data is shown in Fig. |7| 
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Figure 3: The new HERA data on F2 versus x at fixed Q\ 
0.75, g s = 4.019, A = 2. 
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Fig. H-a 



Fig. H-b 



Fig. @-c 



Figure 4: Comparison of solution I (dotted line) and solution II ( solid line) for three values of Q 2 (all other 
information is in the pictures.) 



One can see that with gq which does not increase with energy 
we cannot obtain the good fit of the experimental data for the total 
proton-proton cross section. However, it is clear that the solution 
II gives the description closer to the data in comparison with the 
solution I. The fact that we did not obtain a good description of 
the data for the proton-proton scattering does not look discourag- 
ing to us since we made oversimplified assumption about ^proton- 
colourless dipoles are correct d.o.f. at high energy and Ko(rq) is the 
wave function of the dipole. It should be stressed that the unknown 
mechanism which is different from N=4 SYM and which leads to 
do contribution in Eq. ( 2.17] ) is responsible only for the half of the 
total inelastic cross section at RHIC energies (W = 300 GeV) and 
less than a quarter for the LHC energies. 
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Figure 5: F 2 (x = 0.01; Q 2 ) = 
F™ (Q 2 ) versus Q 2 for solution 
I(dotted line) and solution II(solid line). 
In this figure the errors are shown since 
the values of F™ (Q 2 ) were considered 
as independent fitting parameter. 



4. Conclusions 



As has been mentioned comparing the N=4 SYM prediction with the experimental data we obtain two 
surprising results. First, the N=4 SYM formula gives a good description of the DOIS structure function in 
wide range of Q 2 (Q 2 = 0.85 -f- 60 GeV 2 and x ( x < 0.01). The surprise stems from the fact that Eq. (2.1) 
leads to power-loke dependence on x (F2 oc (l/x) 1_p ) and this power does not depend on Q, Comparing 



10 



4> , Q =0.85 GeV 2 



W=14TeV 



this formula with Fig. one can conclude that dependence of z of Eq. (2.1) as well as on Ins simulates the 
effective power dependance on Q. 

The second surprise is the smallness of g s that generates a very 
small shadowing corrections which we can neglect even at the high- 
est accessible energy:!^ = 14:TeV. p = 0.7 -i- 0.75 means that the 
intercept of the Pomeron is rather small Ajp ~ 0.3 -j- 0.35. We used 
to consider such a small intercept to be typical for the weak cou- 
pling limit (for the BFKL Pomeron). On the other hand, in small 
coupling limit we expect a strong shadowing correction induced by 
the Pomeron interactions ( see Refs. [|ll|, 12, 15, |0]]). Recalling 
that p = 0.7 -r- 0.75 corresponds to A « 7 4- 8 we could expect that 10 
the corrections of the order of 1/A 2 will be small. In this case we 
expect the small shadowing corrections with our fitted small value ° 
of g s . Therefore, we have a dilemma: either the corrections of the 
order of 1/A 2 are large or the shadowing phenomenon is neglidgibly 
small. 




10 



10 



Figure 6: Functions $ 7 (two dashed 



lines), ^proton ( s °lid line 3) and the am- 
plitude of Eq. ( |2.3| ) (solid lines 1 and 2) 
The influence of the corrected b dependence was expected but as functions of z 2 = z for fixed z\ = 
the fact that even with corrected b dependence we have still small lOGeF -1 at two energies W = y/s = 
shadowing corrections was not expected. 20 GeV and W = vs = 14TeV^ Dot- 

ted curve presents Inx4 of Eq. ( 2.16| ) at 
In general we believe this analysis of the experimental data in w_ r — ~\at V 

the framework of N=4 SYM theory gives the useful information on 

the possible scenario what is going on in strong coupling limit at high energy. The picture that arises from 
this analysis is in clear contradiction from the expectation of the high density QCD and because of this it 
could lead to a better understanding the matching between soft (large coupling) and hard (small coupling) 
processes in QCD. 
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Figure 7: The comparison of the total 
cross section for proton-proton high energy 
scattering with N=4 SYM predictions for 
solutions I (dotted line) and II (solid line). 
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